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Abstract: We present the two-loop helicity amplitudes for the scattering of massless quarks and 
massless gauge bosons in QCD. We use projector techniques to compute the coefficients of the general 
tensor describing the two-quark two-boson amplitude in conventional dimensional regularisation and 
use these coefficients to derive the helicity amplitudes in the 't Hooft-Veltman scheme. The structure of 
the infrared divergences agrees with that predicted by Catani while expressions for the finite remainders 
for qq — > gg, qq — ► 57, qq — > 77 and the processes related by crossing symmetry are given in terms 
of logarithms and polylogarithms that are real in the physical region. We have checked that the 
interference of tree and two-loop helicity amplitudes, summed over helicities and colours, reproduces 
the previous results for the finite remainders for interference of tree and two-loop amplitudes given in 
Refs. [1,2]. We also find agreement with the two-loop helicity amplitudes for quark-gluon scattering 
presented in Ref. [3] . 
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1. Introduction 

QCD perturbation theory has proven to be a reliable tool with which to generate theoretical predictions 
of hard scattering processes that can be compared with experimental results from electron-positron 
annihilation, electron-proton and hadron-hadron scattering. In such studies, the theoretical analysis of 
an observable at next-to-leading order (NLO) in a s has become a standard candle, the best theoretical 
prediction to compare against experiment. These NLO calculations are also used to provide an estimate 
of the uncalculated theoretical uncertainties via the variation in the renormalisation (and factorisation) 
scale. 

However, the large amounts of data soon to become available from for example, high-energy 
hadronic collider runs place a demand for more accurate theoretical predictions. To go beyond NLO 
to the next order in QCD perturbation theory, next-to-next-to- leading order (NNLO), has a direct 
impact on several aspects of the theoretical predictions. For example [4], the renormalisation (and 
factorization) scale uncertainty is significantly diminished, there is a better matching of the parton 
level jet algorithm used in theory and the hadron level jet algorithm used by experiment and the 
contributions attributed to non-perturbative power corections may also be reduced. 

There are several ingredients that need to be considered to calculate the NNLO corrections for a 
specific 2^2 scattering process, 

• tree-level 2^4 partonic amplitudes 

• one-loop 2^3 partonic amplitudes, 

• two-loop 2^2 partonic amplitudes. 

The tree-level six-point amplitudes [5, 6, 7, 8, 9, 10, 11, 12] and one-loop five point amplitudes [13, 
14, 15] have been known for some time while there has been recent progress in computing two- 
loop amplitudes. Parton distribution functions from a global fit to available data with observables 
computed consistently at NNLO [4] coupled with a complete knowledge of the three-loop splitting 
functions [16, 17, 18, 19, 20, 21, 22, 23, 24] to generate the evolution of the parton distribution 
functions at this order are also required. 

In this particular paper we are concerned with processes involving two massless quarks and two 
massless gauge bosons, 

qq ->• gg, qq gj, qq 77 

and the processes related by time reversal and crossing symmetry. In previous work [1, 2], we computed 
both the self-interference of the one-loop amplitudes and the interference of the tree and two-loop 
amplitudes for these processes in conventional dimensional regularisation (CDR). In that work, the 
matrix-elements were averaged over all external helicities and summed over colours in the high energy 
limit where all quark masses can be neglected. Here, we extend these calculations to now compute 
the helicity amplitudes for these processes. 

The one-loop helicity amplitudes for parton-parton scattering processes have been known for some 
time [25]. More recently, there have been a few two-loop helicity amplitudes calculated within CDR 
with d = 4 — 2e dimensions: 
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• for the gg — ► 77 [26], 77 — > 77 [27, 28] and 55 — > 5(7 [29, 30] scattering processes, the explicit 
contraction of the amplitude with the external polarization vectors was used. 

• for the for the 77 — > 77 [28] and e + e~ — > 3 jets [31] scattering processes, the individual compo- 
nents of the Lorentz-invariant amplitude decomposition were projected out. 

After the explicit helicity structure is obtained via any of these methods, one needs to tackle a 
large number of two-loop integrals with different topologies and tensor structure. It is now well known 
(see for example [32, 33, 34] and references therein), that these integrals can be reduced [35] to a 
handful of master integrals using integration-by-parts [36, 37] and Lorentz-invariance [38] identities. 
For the process at hand, the relevant master integrals are 10 massless four-point functions with all 
legs on-shell which have been calculated in the past few years using a wide variety of methods [39, 40, 
41, 42, 43, 44, 45, 46, 47, 48]. 

In this paper, we follow a similar appoach to extract the helicity structure of the amplitude as the 
one outlined in Refs. [28, 31]. We start with a general tensor that describes the two quark-two boson 
amplitude and use gauge invariance to find relations among the tensor coefficients. These coefficients 
can be extracted using a set of projectors that are valid at any order in perturbation theory. The 
projectors saturate the tensor structure and produce scalar integrals. The methodology of reducing 
these scalar integrals to master integrals is exactly the same as for the calculation of the spin summed 
matrix elements [1, 2]. Note that the tensor coefficients are evaluated in CDR. 

Then, we provide the perturbative expansion for the one- and two-loop tensor coefficients and 
remove ultraviolet (UV) divergences at each order in a s by renormalisation within the MS scheme. 
The infrared (IR) divergent structure is shown to agree with the one obtained in Refs. [49]. Note that 
in extracting the coefficients of the general tensor structure, conventional dimensional regularisation 
is preserved. 

Finally, the helicity amplitudes are obtained from the general tensor using standard helicity tech- 
niques and within the t'Hooft-Veltman (HV) scheme where the external states are four-dimensional 
and the internal states are kept in 4-2e-dimensions [50, 51]. The finite helicity amplitudes are the main 
new results presented in this paper and we give explicit analytic expressions valid for each helicity 
configuration and process in terms of logarithms and polylogarithms that are real in the physical 
domain. 

During the writing of this paper, we became aware of an independent calculation of the two-loop 
helicity amplitudes for quark-gluon scattering [3]. This paper [3] uses a rather different method to 
extract the helicity amplitudes in a variety of schemes. Where the two calculations overlap (in the 
HV scheme), we have compared the finite remainder terms for the two-loop qq — > gg processes and 
find complete agreement. 

The paper is organised as follows. In Section 2 we introduce notation and discuss the general tensor 
structure of the two quark-two boson amplitude. Gauge invariance is used to reduce the number of 
independent tensor structures. We give projectors to isolate the various coefficients and relate the 
individual coefficients to the interference of amplitudes. Section 3 shows the perturbative expansion 
for the tensor coefficients computed within CDR and presents an analysis of their UV and IR divergent 
structure. The coefficients are renormalised in the MS scheme to remove all UV divergences, while the 
IR poles are predicted using the results of Ref. [49] . Expressions for the one-loop tensor coefficients 
to all orders in e appear in Section 3.2 in terms of one-loop bubble integral and the one-loop box 
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integral in d = 6 — 2e. The relationship between the tensor coefficients and the helicity amplitudes 
in the HV scheme as described in Section 4. The UV and IR structure of the helicity amplitudes is 
exactly the same as for the tensor coefficients, and finite O (e°) remainders are obtained by subtracting 
the predicted pole structure from the explicit calculation both at one-loop and two-loop. The main 
results of this paper are the finite remainders of the two-loop helicity amplitudes. The expressions 
are rather lengthy and we confine ourselves to presenting the amplitudes for the qq — > gg process in 
Appendix A. For completeness, the finite one-loop remainders are collected in Appendix B. Results 
for the processes related by crossing symmetry and those processes involving photons can be obtained 
from the authors. Finally, in Section 6 we present a discussion of our results and conclude. 



2. Notation 



The generic process we consider in this paper is 

q(pi, Ai) + q(p 2 , A 2 ) + b 3 ( P3 , A 3 ) + b 4 (p4, A 4 ), (2.1) 

where the bosons 63 and 64 are either gluons or photons. All particles are light-like with momenta 
satisfying, 

o-p? + p£ + i£ + p£, p1 = o. (2.2) 

Physical processes are obtained by crossing particles into the initial state. The associated Mandelstam 
variables are given by 

■S12 = (Pi +P2) 2 , s 2 3 = (P2 +P3) 2 , si3 = (Pi +P3) 2 , sn + S23 + S13 = 0. (2.3) 

We work in conventional dimensional regularisation treating all external quark and boson states 
in d dimensions and renormalise the ultraviolet divergences in the MS scheme. The bare coupling «o 
is related to the running coupling a s = a s (fi 2 ) at renormalisation scale fi, by 



where 



(!) + (§-§) (S) 2 +^) 



(2.4) 



S e = (47r) e e - £7 , 7 = 0.5772 . . . = Euler constant (2.5) 

is the typical phase-space volume factor in d = 4 — 2e dimensions and /2q is the mass parameter 
introduced in dimensional regularisation [52, 53, 50, 51] to maintain a dimensionless coupling in the 
bare QCD Lagrangian density. 

The first two coefficients of the QCD beta function, /?o and j3\, for Np (massless) quark flavours, 

are 

. 11C A - 4T R N F 17C 2 A -10C A T R N F -6C F T R N F 
Po = q , Pi = , (i.bj 

where iV is the number of colours and 

N 2 -1 1 
Cf = ^T' ° A = N ' TR= r (2 ' 7) 
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2.1 The general tensor 

The most general tensor structure for the amplitude, \M), is 



\M) = u(pi)fou(p 2 ) (Ane3.p1 e4.p1 + A12e3.p1 e 4 .p 2 ) 
+ u(pi)p , 3u(p 2 ) (A 21 e 3 .p 2 e4.p1 + A 22 e3.p2 e4.p1) 
+ u(pi)i4u(p 2 ) (Bxe3.p1 + B 2 e 3 .p 2 ) 

+ u(pi)feu(p 2 ) (C1e4.p1 + C 2 €4.P2) 
+ Diu(pi)^s^4u(p 2 ) 

+ D 2 u(pi)^4hh u (P2) , (2.8) 

where the constraints £3-^3 = 64-^*4 = have been applied. All coefficients depend on the identity of the 
gauge bosons and are functions of si 3 and S23 (and implicitly si 2 = —S13 — S23) where Sij = (pi +pj) 2 . 
This tensor structure is a priori d-dimensional, since the dimensionality of the external states has 
not yet been specified. The amplitude \M) obeys the Ward identity when one of the gauge boson 
polarization vectors is replaced with the corresponding boson momentum, 

\M)(e 3 ^ P 3) = \M)(e4^P4) = 0. (2.9) 

These constraints yield relations amongst the 10 distinct tensor structures and applying these identities 
gives the gauge-invariant form of the tensor, 

\M) = Aii(si 3 ,S23) Tu + A l2 (si3,S 2 3)Ti 2 
+ -42l(si3, S23) T 2 l + A 22 (S13, S23) T 2 2 

+ 5(513, s 23 )T, (2.10) 

where Ajj and B are five gauge-independent functions and the tensor structures T/j and T are given 
by 

TlJ = U(pi)p , 3u(p 2 )e3.pi€4.pj - ^Yu(pi)i3u(p 2 )€4.pj - ^u(pi)^ 3 ^4u(p 2 ), (2.11) 

T 2 j = u(pi)i>3u(p 2 )e3.p 2 e4.pj - -^-u(Pi)fe u (P2)e4-PJ - ^u(p 1 )^ 3 ^ s u(p 2 ), (2.12) 
T = s 2 3 (u(pi)i4u(p 2 )e3.p 1 + ^u(pi)fefofau(p 2 

- S13 (u(pi)i4u(p 2 )e3.p 2 + ^u{jpi)i4hh u (P2)^ , (2-13) 

with S14 = S23- Each of the tensor structures satisfies the Ward identities. The tensor coefficients 
are further related by symmetry under the interchange of the momenta (and colour labels) of the two 
gauge bosons, 

^2l(si3,S23) = -^12(S23,S1 3 ), 
-422(513,523) = -^11(S23,S13), 

B(S 13 ,S 2 3) = B(s 23 , si 3 ). (2.14) 
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2.2 Projectors for the tensor coefficients 

The coefficients Ajj and B may be easily extracted from a Feynman diagram calculation, using 
projectors such that 

£>P0 \M) = X(s 13 ,s 23 ). (2.15) 

spins 

The explicit forms for the projectors in d space-time dimensions are, 

P(Am) = d T f + ^ - ^ T f + 1 r t 

V 11 J 2s 2 2 s 13 s 23 (d-3) 11 2s 2 12 s 13S23 (d- 3) 22 2s 2 2 s 13 s 23 (d - 3) 

(2si3-(ri-4)s 2 3) T t , (2823 4)313) t 

+ 2 S 2 2S 2 3S23(d _ 3) 12+ 2s 2 12 s 13 s 2 3 (d-3) 21 ' 



= (2 8 i3-(d-4)a 23 ) t (2ai 3 -(d- 4)823) t _ (2ai 3 - (d - 4)g 23 ) t 
28 2 2 8 2 3 8 23 (d-3) 11 2 S 2 2 s 2 3 s 23 (d- 3) 22 2 S 2 2 s 2 3 s 23 (d-3)(d-4) 



(4(d - 3)s 2 l3 + 4(d - 4)s 13 s 23 + d(d - 4)s 2 3 ) f 
2sf 2 ^23(d-3)(d-4) 12 

(2(rf - 4)g? 2 + (d - 2)(rf - 6)^13523) T t 

28 2 2 3 2 3 8 2 3 (d-3)(d-4) 



-p(_g) = T^ T^ H 

V 7 2 S f 2Sl3 s 23 (d - 3) 11 2^^3523^- 3) 22 2s 2 2 8i 3 8 23 (d-4) 

(2si 3 -(d- 4)8 23 ) (2823 ~(d- 4)813) 



"J 1 , + 



2s 2 2 8 2 3 s 2 3(d - 3)(d - 4)~ 12 28 2 2 si 3 8 2 3 ((i - 3)(d - 4) 2i ' 

The projectors for A21 and yl 22 are obtained by the permutations 

V(A 22 ) = V(A n ) {si 3 <- 823, rfj <- Ti,rf 2 <- T^Tt - -Tt} , 

p(^2i) = r(A 12 ) [ Sl3 <- S23 ,tI <- ri,rf 2 <- Tj^Tt - -Tt} . 

2.3 Reconstruction of the matrix elements 



(2.16) 



(2.17) 



Starting from the tensor coefficients, the interference (or square) of amplitudes can be easily extracted 
in a variety of schemes. In the HV scheme, the external particles are treated in 4-dimensions and we 
find that, 

3 

(M\M) = (2B - An + Ai 2 ) (2B - A u + A 12 ^ 

spins 



si 3 s 23 sf 2 



2 

(AnA^ + A 22 A\ 2 ) . (2.18) 

Note that the tensor coefficients are in general matrices in colour space and the colour sum has still 
to be performed in evaluating, for example, yliii^. In CDR, additional terms of order e and e 2 are 
straightforwardly generated. The resulting formula is, however, rather lengthy and we do not give it 
here. 
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3. The perturbative expansion of the tensor coefficients 



The three processes that we wish to consider in detail are, 



-> qqgg, 
-> qqgj, 
-> gg77. 



(3.1) 
(3.2) 
(3.3) 



The amplitude for the process P (where T 3 = gg, gj and 77) has the perturbative expansion of 
the form, 



\M 



\M 



■p,(o)\ 



+ 



2^ 



l^ (1) ) + (g) 2 |^' (2) ) + ^) 



where the overall coupling factor is given by 

c v = At: aKael) 1 -* , 



(3.4) 



(3.5) 



with q = 1,1/2,0 for V = 55,(77,77 respectively. \Ai' p ^ n ' > ) denotes the n-loop contribution to the 
amplitude. Denoting the tensor coefficients as X^ (for X = Ajj and B), we find the perturbative 
expansions of the form 



X 



u 



2vr/ 



(3.6) 



where the dependence on the kinematic scales is implicit. The n-loop perturbative coefficients are 
vectors in colour space and can be further decomposed as, 



X99,(n) _ ^2 C 99 '^ X 99 '^'^, 



i=l 



where 
and 



X9~f,(n) — £37,[1] _X"S7>(»»),[1] 
£77, [1] _ g Jf 77 '^'W, 

£33, [1] _ ^ a3 T a4 )- • , £99, [ 2 ] = (y»4y a 3^ ^ £33, [3] _ £ a 3a4J 



£37,[1] — (T" a 3). . ; £77, [1] — § 



(3.7) 

(3.8) 
(3.9) 



Index i denotes a colour label in the fundamental representation while index a is in the adjoint. Note 
that the normalisation of the colour matrices is given by Tr(T ai T a2 ) = ^S aia ' 2 . 
Symmetry under interchange of the two gluons gives the additional relations, 

A99,(n),[2] i \ _ A gg,{n),[l} ( x 



A 



33,0), [3] . 



J J \S13,S23) = -Ajj [S 2 3,S13) 

B 33,(n),[2] (si3)S23) = B33,W,[1] (S23)S13)) 

5 33,H,[3]( Sl3)S23 ) = B "^ 3 \ S23 , Sl3 ). 



33,0), [3], 



(3.10) 
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For the processes involving one or more photons, 

Ar/ n) ' [1 \ S is,s 23 ) = -A^ n) ' [1 \s 23 ,s 13 )^ 

B^WM(s 13 ,s 23 ) = B^ (n ^(s 23 ,s l3 ), 

£77,(n),[l] (si3)S23) = B77,(n),[l]( S23)Sl3) . (3J1) 

As in previous work, we use QGRAF [54] to generate the Feynman diagrams and then use 
MAPLE [55] and FORM [56] to manipulate the Dirac structures and traceology and then reduce 
the scalar loop integrals to master integrals. 

3.1 Tree level results 

At tree level, none of the Tu structures are present so that, 

A V I j° m {s l3 ,s 23 )=^ (3.12) 
while only the first two colour structures contribute to B"'(°\ 

£^(0),[i]( Sl3)S23 ) = 2_ 

Sl 2 S 23 

B^M(s 13 ,s 23 ) = — , 

S12S13 

£^(0),[3] (si3)S23) = Q (313) 

For the other processes, 

#n>(0),[i]( Sl3) S2a) = i? 77,(o),[i] (si3) S2s) = B ra,(o),[i] (si3) S23 ) + B M>WM(s 13 , s 23 ) 

2 (3.14) 



Sl 3 S 23 



3.2 One-loop results 



At one-loop, all of the tensor and colour structures contribute. For the two-gluon process, we find the 
following (unrenormalised) structure, 

^?f (1),[1] (si3, s 23 ) = C F A(s l3 , s 23 ) + C A9l (s 13 , s 23 ) + (2N F - C A (d - 2)) h(s 12 ), 

A gg,(l),[2], v _ A gg,(l),[l], , 

A ll {Sl 3 ,S 23 ) — -A u [S 23 ,Si 3 ), 

/ <9 , S,(1)>[3] / \ — r\ 

A u v s i3, S23J — U, 

A 9 l 9 2 ' mi] (s 13 ,s 23 ) = C F f 2 (s 13 ,s 23 ) + C A g 2 (s 13 ,s 23 ) + (2N F - C A (d - 2)) h(s 12 ), 

Af 2 ,(1) ' [2] (si3, s 23 ) = C F f 3 (s 13 , s 23 ) + C A g 3 (s 13 , s 23 ) - (2N F - C A {d - 2)) h(s 12 ), 

^• (1M3] ( S 13, S2 3) = h^SM), 

B^ (1 ^( Sl3 ,s 23 ) = C F fr(s l3 ,s 23 ) + C A g A ( Sl3 ,s 23 ), 

B^,(l),[2] (si3)S23) = i^W-M^^), 

BM.W.[3]( ai3j S23 ) = /l2 ( Sl3; S23 ) j (3.15) 
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while for the gluon-photon and photon-photon processes we have, 

^' (1) ' [11 ( Sl3 ,s 23 ) = C F (f 1 (s l3 ,s 23 )-f 1 ( S23 , Sl3 )) 

C 'a 

+ ~2~ (fl(^13, S23) - gi(s23, S13) ~ 4/l 3 (si 3 , S23)) , 

^• (1) ' [11 ( Sl3 , S2 3) = C F (f 2 (s 13 ,s 23 )+ f 3 (s 13 ,s 23 )) 

Ca 

+ — (#2(513, s 23 ) + g 3 {s\ 3 , s 23 ) - 2hi(s 13 , s 23 )) , 
B^W>W(s 13 ,s 23 ) = C F (f 4 (s 13 ,s 23 )+ f 4 (s 23 ,s 13 )) 

C A 

+ ~2~ (54(513, s 23 ) + 54(^23, s 13 ) - 2h 2 (s 13 , s 23 ) - 2h 3 (s 13 , s 23 )) 



and, 



(3.16) 



AlT mi \s 13 ,s 23 ) = C F (fi(s 13 ,s 23 )-f 1 (s 23 ,s 13 )), 
A^ {lUl] (s 13 ,s 23 ) = C F (f 2 (s 13 ,s 23 )+ f 3 (s 13 ,s 23 )), 

B^M(S 13 ,S 23 ) = C F (f A (s l3 ,S 23 ) + / 4 (S23,*13)) • (3.17) 

The functions fi, gi, hi and h can be written in terms of the one-loop box integral in d = 6 — 2e 
dimensions, Box 6 (sjj, s^), and the one-loop bubble, Bub(sy). For all orders in e = (4 — d)/2, we find 
that, 

2(d-4) 
5? 2 (d-l)(d-2) 



M«i2) = 2 o> Bub( Sl2 ), (3.18) 



while, 



(<*-2)(d-4)„ fi/ , 2(d-4)„,, , 2(d-4)„,, , 
/1 Si 3 ,S23 = ^ ^ Box 6 si2,s 2 3 - - - Bub s 23 " - " Bub (s 12 , 

«13 Si 3 S 23 S 12 S 13 

s 2 23 d(d - 2) - s 13 s 23 (d -2){d- 14) - 2sUd - 10) 6 

/2(Si 3 ,S 23 J = o DOX (512,^25) 

S1 2 S| 3 

2{d - 12)s 13 - 2s 23 d 2s 23 d-4{d-8)s l3 
H o Bub(s 23 ) H 2 Bub(si 2 ), 

Sl2S{ 3 S 12S{ 3 

{d-A){d -10)^3 + (3d 2 -36d + 92) Sl3S2 3-2(d-2) S f 3 6 

M s 13, S23j = DOX (Sl2,Sl3j 

Sl2-Sl 3 S 23 

4si3-4(d-7)s23 v , 2(d-10)s 23 -4si3 D , 
H Bub(si 3 ) H Bub(si 2 ), 

Sl 2 Sl 3 S 23 Sl 2 Sl 3 S 23 

4(d - 3) g f 3 - (d 2 - VSd + 38)g 23 s 13 + (d - 2)s 2 23 6 

/4(S13,S23) = H ~ BoX D (si2,S 23 ) 

Sl 2 Sl 3 S 23 

, 2(d-7)ai3-2a 23 p u , , 2{d - A)s 23 + 8(d - 3)s 13 ^ 
H Bub(s 23 ) H 7^ — 7^ Bub(si 2 ), 

S12S13S23 (d - 4)si 2 Si 3 S 23 

g\{si 3 ,s 23 ) = -/i(si 3 ,s 23 ), 

S23(d-2)(4si 3 + ds 23 ) 6 . . 2s 23 d + 4si 3 
52(si 3 ,s 23 ) = 2 Box ( s i2,«23)H o Bub(s 23 ) 

Sl2Si 3 Sl2Si 3 

2(d- 6)^13 -2dg 23 p , , , 
H 2 Bub(si 2 ), 

S1 2 S| 3 
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53(si3,S23j = Box (si2, S13) H Bub(si 2 ) 

S12S23 -512S23 
2(d-4)s23-4si3 

H Bub(si 3 ), 

•S12S13S23 

54(Sl3,S23j = o B0X D (S12,S2 3 ) 

^•Si2Sl3 

(rf 2 -16d + 40)£i3 -2(d-4) g23 p , , . 2 
— — Bub s 2 3 Bub(si2 , 

(d - 4)si 2 Si 3 S23 Sl2«13 

, , , (d -10)((d- 2)^3 + 2si 3 ) p 6 , , 

«l(si3,S23j = ~ BoX°(si2,S 2 3) 

^•Sl2«13 

(d-4)(d-10)(2 Sl3 + S2 3) p .6, , , 

2si 2 Si3 

(d-10)(2si 3 -(rf-4)s 23 ) p 6 , . 
o Box D (si3, s 23 ), 

2S12S13 

h , , (8(d-3)s 13 -(d-4)(d-10) S23 ) ( T) 6 6 ' 

ra 2 (si3, S23) = 4si 2 s 2 3 I Box°(si2, S23) + Box u (si 3 , S23) 

(8(d _ 3)g23 -(d-4)(d-10)s 13 ) f p 6 6 ' 

H 4si 2 si 3 I Box°(si2,si 3 ) + Box u (si 3 ,S23j 

4 + 3(d-4)2 

; -B0X 6 ( Sl3 ,S23) 



4S12 

4(d-3) / 1 



Bub(s 23 ) + Bub(si 3 ) + Bub(si 2 ) , 



(d-4) V s 12«13 S12S23 ' S23S13 

^3(S13,«23) = 4 /l(si2,S 2 3)- (3.19) 

3.3 Ultraviolet renormalisation 

The renormalisation of the matrix element is carried out by replacing the bare coupling «o with the 
renormalised coupling a s = a s (fi 2 ), evaluated at the renormalisation scale /j, 2 according to Eq. (2.4). 

We denote the i-loop contribution to the unrenormalised tensor coefficients for process V by 
4j W ' un and B^.W,™ 

using the same normalization as for the decomposition of the tensor coefficients 
(3.4); the dependence on (si3,s 2 3) is always understood implicitly. The renormalised coefficients are 
then obtained as 

Ajj - U, 

aT,W - q -l A V,(l),un 
A IJ — D e A IJ ' 

A Vj 2 ) = S- 2 A^} {2),mi - ( g + g-i A^'™ 1 , (3.20) 

and 

B V,(0) = B V,(0),m^ 

B V,(l) = s -l B V,{l),un _ QPo B P,(0),nn^ 

B V,(2) = s -2 B V,{2),nn _ (Q + 1)A) g-l g p,(l),un _lf^_ (<l+ 1 )f 3 o \ B 7>,(0),un 

6 2 V 6 6 / 
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(3.21) 

where q is the number of powers of the strong coupling for process V as described in Eq. (3.5). 
3.4 Infrared behaviour of the tensor coefficients 

After performing ultraviolet renormalisation, the amplitudes still contain singularities which are of 
infrared origin and will be analytically cancelled by those occurring in radiative processes of the 
same order. Catani [49] has shown how to organize the infrared pole structure of the one- and two- 
loop contributions renormalised in the MS scheme in terms of the tree and renormalised one-loop 
amplitudes. The same procedure applies to the tensor coefficients. 

In particular, the infrared behaviour of the one-loop coefficients is given by 

A V,(1) _ ,7>,(l),finite 

B V,(1) = jr7>,(l)( e ) B 7>,(0) + B P,(l),finite ) ( 3 _ 22) 

while the two-loop singularity structure is 

A Vf) = J^(l) (e) ^ ; « + ^), fi nite ) 

B V,(2) = jV,(2) (e)jB P,(0) + jV,(l) (e)j g(l) + B (2),finite j (3 _ 23) 
where ^'^(e) is given by, 

j^,(2) (e) = _ 1^,(D (e )I^« (e ) - ^.(i) (e) + e -er EgLz^ ( ^ + k) I v ^{2e) 

— C J. I J. € j \ £ / 

+H v + 2 \e) (3.24) 

and the constant K is 

«=(%-4)ca- ^TrNf. (3.25) 



^18 6 J 9 
In QCD the /^'^(e) operator is given by [49] 



^' (1) W = ^r^E^)E T *- T i( 

i j^i 



2pi • pj 



(3.26) 



where the sum runs over the pairs of external coloured particles involved in process V and the singular 
function is 

zA m9 (e) = - + 7i -, 7(M - = -, lg = ^- (3.27) 
The colour charge algebra is defined as 



Tj-Ti if i ^ j, 
T 2 = d otherwise, 



with the Casimir C\ = Cp (Cj = Ca) if parton i is a quark (gluon) and where the colour charge Tj is 
^cb (~^cb) f° r a ci uar k (anti-quark) and if ca b for a gluon. 
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Then, given the colour basis for the two-quark two-gluon process presented in Eq. (3.7), we have 
the following infrared-singularity operator 



JOT.(l)( e ) 



=£7 



T(l - e) 



f A(e,s 12 ,s 13 ,S2 3 ) 

A(e, si2, «23, S13) 

y B(e, S12, S13, S23) B(e, s 12 , S23, S13) 



B(e,S 2 3,Stf,Sl2 
B{e,Si3,S 2 3,Sl2 

C(e,s 12 ) 



\ 



where 



A(e, S12, S23, S13) = - 



B(e, S12, S23, S13) 
C(e,s 12 ) 

Similarly, for the other processes we have 

JS7,(l)( e ) 





2r(l-e) 



-iV 



1 



3 



2Ne 



S\2 
2 



.E. 
S23 



S13 



and 



J77,(l)( e ) 




(3.28) 



(3.29) 
(3.30) 
(3.31) 



(3.32) 



(3.33) 



In Eq. (3.23) the function H^'^ contains poles of 0(1/ e) and is process and renormalisation 
scheme dependent. The specific structure of this term was not given in Ref. [49]. However, there it 
was pointed out that in the case of a QCD process that involves a qq pair, the function H^'^ is 
related to the quark electromagnetic form factor. Now, we know that the form factor is not merely 
related, but is the backbone of the singular structure of multi-loop QCD amplitudes [57]. 

It was only after the direct Feynman diagram evaluation of two-loop matrix elements [58, 59, 60, 
1, 61, 2], that the colour uncorrelated part of the H^'^ function for all 2 — > 2 partonic processes was 
fully identified. More precisely, it was found that each external coloured leg in the partonic process 
contributes independently to the matrix element given by 



{M \H V ^\M ) 



on 



4eT(l 



-# p <( 2 )(.Mo|A"*o> 



(3.34) 
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with 



H v,(2 ) = nqH (2 ) + ngH m (335) 

where n g is the number of external gluons, n q is the number of external quarks and anti-quarks for 
the process V and 

+ (-U + ^) TrNfCf ' (3 ' 36) 

Cs , 5 , 11 _ 2 \ ^ 2 



+ ( v Y + 12 + 144 7r J C - < 3 - 37 > 

Already in Ref. [49], Catani pointed out the possibility of H^'^ having additional non-trivial 
colour correlations which could be explored by calculating the hard scattering amplitude, rather than 
the interference of two-loop and tree amplitudes. 1 Such colour correlations were found for the four 
gluon interaction, which is the only partonic process (involving four coloured particles) for which 
the two-loop helicity amplitudes [30] exist. According to Ref. [30], the colour correlations form a 
non-trivial matrix in colour space and is the same for all helicities. 

For the processes we study, we find that the colour correlations are indeed helicity independent 
and are given by 

^ <2, (^4lf^)(~) 2 '(^• ,2)l + fiP • ,2, ) < 3 - 38 ' 
with H V '( 2 } as given in Eq. (3.35) and 

H99,(2) = _ 4 i f abc T a[i] T m T cm ln ( zM ln ( Z^l) i n ( Zfll) j (3.39) 

\-S23j \-Sl3j \-Sl2j 

H S7 '( 2 ) = 0, (3.40) 
H 77 '( 2 ) = 0. (3.41) 

These results are in complete agreement with the results found in Ref. [30] and show that the colour 
correlations are also independent of the external particles. This also agrees with the analysis of 
Ref. [57] where it was shown that the contributions to H^'( 2 ) are of soft nature and directly related 
to soft anomalous dimension matrices obtained by studying the evolution of colour exchange in QCD 
hard scattering (see [62] and references therein). It is worth mentioning that these type of colour 
correlations between three different partons also appeared in the analysis of higher-order contributions 
to the soft-gluon current in Ref. [63] . 



lr The sum over colours causes the colour correlations to vanish 
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Note that the matrix f abc T a ^T b ^T c ^ collects all colour weights that result from the analysis 
of triple gluon vertices connecting any three different external partons, projected into the colour basis 
given in Eq. (3.7), so that for the two quark-two gluon process, 



i f abc T a ^T h ^T^ = -- 



( 2 2 4N ^ 

-2 -2 -AN 
y-N N j 



(3.42) 



Given the tree level results presented in Section 3.1 and the divergent structure presented in Eq.(3.23), 
we can see that in fact the third column of this matrix is not needed at O (1/e), since the S^^di^ 
colour structure is not present at tree level. 

4. Helicity amplitudes 

The helicity amplitudes |-M)a 1 a 2 a 3 a 4 can be obtained from the general d-dimensional tensor of Eq. (2.10) 
using standard four-dimensional helicity techniques corresponding to working in the 't Hooft-Veltman 
scheme. Conventionally, the two helicity states of a four-dimensional light-like spinor ip(p) are denoted 
by, 

^±(p) = ^(l±75Mp), (4- 1 ) 

with the further notation, 



\p±) = il>±(p), (P±|=^±(p). (4.2) 

Particles may thus be crossed to the initial state by reversing the sign of the helicity. The basic 
quantity is the spinor product, 

(pq) = {p~\q+), [pq] = [p + IH> ( 4 - 3 ) 

such that 

(pq)[qp] = 2p.q. (4.4) 

In this notation, the polarization vector of the out-going light-like gauge boson of momentum p is 
defined with respect to a reference momentum k as, 

eHp,k) = ± iP / - ±h ^ ±) . (4.5) 

Fixing the helicity of the (final-state) quark to be positive, Ai = + we find that the helicity 
amplitudes, upto a helicity dependent phase, are given by 



I ( i\ v S 12 S 13 S 23 , , \ 
OC -5! Aii(si3,a 2 3), 



r-3 

\M)+-+- oc X^23 (2jB(si3; S23) + Al2 (s 13 , s 23 ) - A u (s 13 , s 23 )) , 
\M)+—+ oc ^^23 (2B( Sl3 , s 23 ) + A 22 (s 13 , s 23 ) - A 21 (s 13 , s 23 )) , 
\M)+— oc - ^ 13823 A 22 («i3, -23). (4-6) 
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By inspection of Eqs. (3.10) and (3.11) we see that, as expected, there are two independent helicity 
amplitudes and that the other two helicity amplitudes can be obtained by Bose symmetry and, 

\M)+—+ = \M)+-+-(s 13 s 23 ), 

\M)+ = \M)+- ++ (s 13 <- 823). (4.7) 

Helicity amplitudes where the quark helicity is negative are obtained through parity, 

|A<)_ + _ A3 _ A4 = (|.M) + _ A 3 A4 )*. (4.8) 

By inspection, we see that "squaring" and summing over the helicities generates the interference of 
amplitudes given in Eq. (2.18). Thus, the helicity amplitudes can be obtained directly from a Feynman 
diagram calculation by projecting out the individual tensor coefficients and taking the appropriate 
linear combination. 

Just as for the individual tensor coefficients, the helicity amplitudes for each process V can be 
perturbatively decomposed as, 



\M V ) 



A1A2A3A4 

2 



|-M P ' (0) )a 1 a 2 a 3 a 4 + g) \M P ^) XM + (g) |M^ 2 )> Al W4 + 0{al) 

(4.9) 

where the overall coupling factor is given by Eq. (3.5). As in Eqs. (3.20) and (3.21), the renormalised 
helicity amplitudes are obtained from the unrenormalised amplitudes by 

|M 7 '' (0) )a 1 a 2 A3A 4 = IM^'^Kw,, 

\M P ^) XMiM = 57 1 |-M P ' (1) ' Un )A 1 A 2 A 3 A 4 " ^|^°)'-> AlA2A3A4 , 

\M )A 1 A 2 A 3 A 4 , 



e 2 



(4.10) 



while the infrared singularity structure is given by, 

|X P ' (1) )A 1 A 2 AaA 4 = ^' (1) ( e )|M P ' (0) )A 1 A 2 A 3 A 4 + l-M^'W^Ax^, (4.11) 
l-M^A.w, = I V ^\e)\M^) XlMXzM +I V ^{e)\M p ^) XlX2X3Xi 

+|M P ' (2) ' finitc >A 1 A 2 A 3 A 4 , (4.12) 

where the operators J^'W are the process dependent matrices in colour space given in Section 3.4. 

Similarly, the n-loop helicity amplitudes are vectors in colour space and can be further decomposed 
as in Eq. (3.7) 

3 

\M^) XlX2MXi = £ C " M \M^ n ^) XlX2XsX4 , 
i=i 

\ M rr>W) = C^'W |^ 37 ' {ri) ' [11 ) AlA2 A 3 A 4 , 

\ M ^,(n) } = C 77,[i] |Mrr,(n),[i] )AiAaAsA4 . (4 _ 13) 
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4.1 Helicity amplitudes for physical processes 

The physically relevant scattering amplitudes we wish to describe are the quark-antiquark annihilation 
processes, 



?(P2, +) + q(pi, ") b 3 {p 3 , +) + 6 4 (P4, +), (4.14) 

?(P2, +) + ?(pi, ") &s(P3, +) + & 4 (P4, -), (4.15) 

?(P2, +) + <l(pi, ") h( P3 , -) + fe 4 (P4, +), (4.16) 

?(P2, +) + g(pi, -) &3(P3, ") + &4(P4, -), (4.17) 

and the quark-boson scattering processes, 

q(p2, +) + 6 3 (pi, +) -► g(p3, +) + & 4 (P4, +), (4.18) 

?(P2, +) + 6 3 (P1, +) 9(P3, +) + MP4, -), (4.19) 

q(P2, +) + 6 3 (P1, ") ?(P3, +) + &4(P4, +), (4.20) 

?(P2, +) + 6 3 (P1, ") ~> ?(P3, +) + &4(P4, "), (4.21) 

?(P2, +) + 6 4 (pi, +) & 3 (P3, +) + 9(P4, +), (4.22) 

?(P2, +) + &4(Pl, +) & 3 (P3, ") + 9(P4, +), (4-23) 

q(p2, +) + 6 4 (pi, -) 6 3 (P3, +) + 9(P4, +), (4.24) 

g(P2, +) + &4(pi, -) & 3 (p 3 , -) + g(P4, +)• (4.25) 



where the bosons 63 or 64 are gluons or photons. Processes where the initial state quark has negative 
helicity are obtained by a parity transformation, while processes with two bosons in the initial state 
are given by time reversal. There are additional symmetries when the gauge bosons are identical and 
further symmetries between the various colour components for the two gluon process. 

For convenience, we denote the processes of Eqs. (4.14)— (4.17) as the s-channel, the processes of 
Eqs. (4.18)— (4.21) as the u-channel and the processes of Eqs. (4.22)— (4.25) as the t-channel. s, t and 
u are the usual Mandelstam variables, s = s\2 > 0, t = S23 < 0, u = S13 < and s + 1 + u = 0. 

To present the helicity amplitudes for the various processes, it is convenient to organise the 
amplitude where particles 1 and 2 are in the initial state in terms of a spinor factor S and the colour 
factor C so that, 

|^' (n) )A 2 A 1 A 3 A 4 = £ sSIa 3 A 4 >< Cc M X I^T' (n) ' W )A 2 A 1 A 3 A 4 , (4-26) 

i 

where c (= s, t, u) denotes the channel and the sum runs over the number of colour structures: three 
for the two-gluon processes and one otherwise. 
Explicitly, we have 

Q99,{1] — ^rpa 3 rpa 4 y 
Q99,[2] _ ^rpa 4 rpa s y ^ 

C gg,n = 8 aiai 5 ili2J (4.27) 
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while for the processes with photons we have, 



Cf'W = (T-) ilia , 

Cr>W = S ili2 , (4.28) 

The u- and i-channel colour factors are obtained by crossing symmetry, 

Cl^=C v s ^ (l<-3), C V t ^=C v s ^ (l<-4). (4.29) 

There is considerable freedom in the definition of the helicity dependent spinor factors. For the 
s-channel, for process (4.15) we choose, 

[!] (24) 3 (14) 
°*+-+- (12) (23) (34) (41)' 

e M ,[2] ■ (24) 3 (14) 

s+ - + - (12)(24)(43)(31)' 



c99.[3] _ (TS7 _ C77 _ . (24) 3 (14) 

5 s+ _ + „ - - - -* (13)(14)(23)(24) , ( '-<■..' 



while for process (4.14) we choose, 



P , [{] _ (23) [13] 



The factors for the other s-channel processes, (4.16,4.17) are obtained by Bose symmetry (exchanging 
labels 3 and 4). Note that for the two-gluon process, this symmetry also exchanges the first two colour 
structures, [1] and [2]. The spinor prefactors for the u- and t— channels are also obtained by crossing 
symmetry, 

5 Sa 3 A 1 A 4 = €£iiXsX 4 (Pi ~ P3)' ^SIasAx = < S 21 1 A 3 A 4 ~ ^ ( 4 - 32 ) 

4.2 Tree level helicity amplitudes 

At tree-level, inserting the results for the tensor coefficients for the various processes given in Sec- 
tion 3.1, the only non-vanishing helicity amplitudes in the s-channel are, 

ffigg ,(0),[l]) + _ + _ = .(0),[l]) + __ + = _ 2) 

|jC?f >(°)'[ 2 ]) + _ + _ = \Mf •(°)'[ 2 ]) + __+ = -2, 

'(°))+_+_ = |.M^' (0) > + + = 2, 

l-Mr ,(0) )+-+- = |A? 77 ' (0) )+— + = 2. (4.33) 

The appropriate spinor factors are obtained from Eq. (4.30) while the colour factors are given in 
Eqs. (4.27) and (4.28). Tree helicity amplitudes for the other channels are obtained by the appropriate 
crossing symmetry. 
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4.3 One-loop helicity amplitudes 

The one-loop helicity amplitudes are straightforwardly obtained by inserting the all-orders results for 
the tensor coefficients given in Section 3.2 into Eq. (4.6), expanding the one-loop bubble and box 
integrals around e = 0, and renormalising according to Eq. (4.10). Explicit expansions in e for the 
d = 4 — 2e one-loop bubble and d = 6 — 2e one-loop box graphs can be found, for example, in the 
Appendix of Ref. [1]. 

The finite remainders of the one-loop amplitudes defined through Eq. (4.11) for the two-gluon 
process can be decomposed according to the number of colours and massless quark flavours, 



\Jy[99,(X), [i],fi™tej _ ^433,(1), M _|_ _L_g»S,(l),[»] + N F (J99,W,[<] _ fog MM ^ 



(4.34) 



for i = 1 , 2 and 

^J^(99,W,[3], finite^ _ jjgg,(l),[3] _|_ jggg,(l),[3] (4.35) 



N 

for i = 3 and where 



and /3o is defined in Eq. (2.6). For the other two processes we have, 

\M9lM, [l],finitej _ j^g-y M ,W + — BgT'W'M + NpC^^^ - -foS \M g c y '^^}, 



N 



1^77,(1), [I], finite^ _ ]y j^rY,(l),[l] + ±_ ^77,(1), [1] _)_ ^(777,(1), [1] _ 



(4.37) 



(4.38) 



For clarity, the dependence on the helicities has been suppressed. Explicit expressions for A, . . . , E for 
qq — ► 55 in the physical region, s = S12 > 0, t = S23 < and n = S13 < 0, are given in Appendix B. 
Explicit results for qq — > (75 scattering in the physical region, s = S12 > 0, t = S23 < and u = S13 < 
are given in Appendix B. Formulae for the other channels and processes are available from the authors 
as FORM outputs. We have checked that the finite one- loop helicity amplitudes presented here are in 
agreement with the helicity amplitudes for qq — > gg (and crossed processes) given for the HV scheme 
in Ref. [3]. 

4.4 Two-loop helicity amplitudes 

The main results of this paper are the two-loop amplitudes for processes (4.14)-(4.25). As in the 
tree and one-loop cases, the helicity amplitudes can be directly extracted from the appropriate linear 
combination, Eq. (4.6), of unrenormalised two- loop tensor coefficients obtained by direct evaluation 
of the projectors given in Section 2.2 acting on two-loop graphs. Renormalisation is achieved via 
Eq. (4.10) and the two-loop master integrals expanded around e = 0. The finite remainder of the two- 
loop amplitudes are defined through Eq. (4.12) and, for each process can be decomposed according to 
the number of colours and massless quark flavours, 

\M99M, [i],finitej _ jy2 j^gg ,(2) ,[i] _|_ _gss,(2),[i] + J_ 2 + NN F £)99M,\i} 
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+ A ^K99,(2),[i] , N 2 F gg,(2),[i] , ^ G ra,(2) ) [i] 

-2(5 S \M9 9 ,m,[i\Jinite^ _ + ffi^ |^(OT,(0),[t]^ (4,39) 

for i = 1 , 2 and 

|j^ff,(2),[3],/irote) = NH 99,(2),[3] + }_jgg,(2),[3] + Np jgg,(2),[3] 

+ ^L K gg,(2),m + ^1 rss>(2),[3] 
^7V 2 c ^ JV c 



-2/3 S \M° 9 ' W ' [3] ' finite ) - (0iS + PiS 2 ) \M 9 C 9 ' {0) ' [3] ). 

For the other processes we have, 

^gj,(2),[l],finitej = ^2^97,(2), [1] + £37,(2), [1] + _L q?7,(2),[l] + NNpD 91 ^' [1] 



(4.40) 



■ 3 /3o^|-Mr' (1) ' [1] ' /inite ) - f 3 /?o^ 2 ) |<Mf" (0) ' [1I >> ( I. ID 

2 V 2 8 



1^77,(2), [l],finitej = ^2^77,(2), [1] + £77,(2),[1] + _L £77,(2), [1] + ATAr F £)77,(2),[l] 

+ ^£Tr,(2),[i] + ^Q^c F ^ 2 )'W -/3bS|^>«>W> /<nite ), 



(4.42) 



where Q g is the electric charge of a massless quark of flavour q. 

As in the previous section, the dependence on the helicities has been suppressed. Explicit expres- 
sions for A, . . . ,L for qq — > process in the physical region, s = S12 > 0, t = S23 < and u = S13 < 0, 
are given in Appendix A. Formulae for the other channels and processes are available from the authors 
as FORM outputs. We have checked that the finite two-loop helicity amplitudes presented here are in 
agreement with the helicity amplitudes for qq — > gg (and crossed processes) given for the HV scheme 
in Ref. [3]. 

5. Comparison with previous results 

In recovering the square, or interference, of amplitudes in a particular channel c, the contraction of 
the n-loop colour vector with a conjugate m-loop colour vector (M P ' <ym ' > \ obeys the rule 

(M r ' {m) \M r ' ( - n) ) 

= E E \^' (n %x 2 x s xJM^) XlX2X3Xi (5.1) 

helicities colours 

where |-M c : ''' n ^)AiA2A3A4 contains all of the spinor and colour information for a particular helicity. 
Expanding the colour and spinor factors C and S, and dropping the explicit dependence on the 
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helicities, we see that, 

(M v ' {m) \M p ' {n) ) 

= E E Y. c ? M * c ^ b]s c' [{] * s c' m \M v c ' {nm r\M r c ' {n) ^ 

helicities colours i,j 

= E E^^i^ ,(n) ' [ll )*i^' (n) ' b1 )- 

helicities i,j 

The symmetric colour sum matrix 0C V = ^colours ^ ^c''^ i s given by 



/ jv 2 -i _j_ 1 \ 



2Af 2W 

1 JV 2 -1 

2iV 2N 

1 1 2iV 



1 



(5.2) 



(5.3) 



and, 



CC 91 = 



N 1 



OCT 1 = N. 



(5.4) 



In the s-channel, the process and helicity dependent spinor matrices SS^ = S^'^* S^'^ are given by, 



u 



u 2 ut us 
tu t 2 ts 
su ts s 2 



•^3 + - + - 



CC77 



U 



(5.5) 



and 



iit 

cc99 _ _ 

°°s+-++ - s 2 



1 1 1 



cc91 



lit 

cell _ _ 



(5.6) 



V 111 / 

Similar matrices for the other helicities are obtained by Bose symmetry, while corresponding matrices 
for the u and f-channels are obtained by crossing symmetry, 



SS^ = -SS^(s <- t), SSf = -SS^(s <- u). 



(5.7) 



The explicit minus sign accounts for the exchange of a single fermion between the initial and final 
states. 

Eqs. (5.1) and (5.2) can be used to recreate the full interference of tree and two-loop amplitudes. 
However, we notice that we can immediately separate the singularities from the finite parts using 
Eq. (4.12). The finite remainder, denoted by J- 2 ^ e (s,t,u) in Refs. [1, 2] is thus given by, 



2Re E E \M r c ^ m T IMP^MJinite^ 

helicities i,j 



(5.8) 



Note that by using helicities, we are implicitly treating the external particle states in 4-dimensions 
- the HV scheme. However, the singular contributions to the renormalised amplitudes defined in 
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Eqs. (4.11), (4.12) are given relative to the tree amplitude and therefore any additional terms that 
would be produced in CDR are automatically removed. Treating the external states differently for 
both the n-loop and tree amplitudes does not alter the finite contribution. Direct application of 
Eq. (5.1) for the real part of the tree and two-loop interference (n = and m = 2) in the HV scheme 
therefore generates the same finite parts given in section 3.2 of Ref. [1] for quark-gluon scattering and 
in Section 4 and Appendix A of Ref. [2] obtained using CDR. 

Similarly, the square of one-loop graphs is obtained using n = 1 and m = 1 in Eq. (5.1). In 
principle, the one-loop amplitude should be expanded through to 0(e 2 ). However, as noted in [1, 
30, 64], the finite self interference only contains logarithms (up to the fourth power), but does not 
contain triple and quartic poly logs. These terms naturally arise when the (9(e) and 0(e 2 ) terms 
in the expansion of the one-loop box are multiplied by singularities from the conjugate amplitude. 
Such terms appear in the Catani pole structure, where l-M*- 1 -*) does need to be 

expanded through to 0(e 2 ). The finite parts of the renormalised one-loop self interference remaining 
after Catani's prediction for the singularities has been removed should therefore be identical to those 
obtained by summing the finite parts using Eq. (5.1). In this instance, the finite remainder denoted 
by ^lnUe^ s i^i u ) m i s obtained from 

^2 X % \M£'W' [i] ' finite )* \M^>W^f inite ). (5.9) 

helicities i,j 

We have checked that this procedure correctly reproduces the coefficients given in Section 5 and 
Appendix B of Ref. [2] for the processes involving one or more photons. Unfortunately, the results of 
Ref. [1] are not presented in quite the same way and a direct comparison is more difficult. 2 However, 
we have checked the finite remainders after the Catani poles have been subtracted, both directly from 
the square of the one-loop amplitude (the CDR scheme) and the one-loop square reconstructed from 
the helicity amplitudes (the HV scheme) with the HV results of Ref. [3] and again find agreement. 

6. Conclusions 

In this paper we have presented the analytic expressions for the one- and two-loop virtual helicity 
amplitudes for quark-gauge boson scattering, qq — ► gg, qq — > 57 and qq — > 77 and those processes 
related by crossing symmetry. These amplitudes were obtained by the construction of cZ-dimensional 
projection operators that extract the coefficients of the most general tensor representation for the 
amplitude, order by order in perturbation theory. Once the renormalised tensor coefficients are known, 
the renormalised helicity amplitudes can be constructed straightforwardly as a linear combination of 
the coefficients. 

At the level of the diagrammatic calculation, we applied the projection operators directly to the 
Feynman diagrams that contribute at each order in a s . The projectors saturate the tensor structure 
of the Feynman diagram and yield a set of scalar integrals that can be related to a basis set of master 
integrals with the application of widely used reduction algorithms. In fact, we are using exactly 
the same tools as we did in [1, 2], except in that case, the projector was merely the conjugate tree 
amplitude. Because the projectors exist in d-dimensions, conventional dimensional regularisation is 

2 In fact there is a misprint in the renormalisation term, TZa, of Eq. 4.6 of Ref. [1]. 
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preserved and there is no ambiguity in dealing with 75. Once the general tensor coefficients are 
determined, the external particles can be treated in 4-dimensions (the 't Hooft-Veltman scheme) and 
standard helicity methods used to construct the helicity amplitudes. 

By summing over helicities and colours, the full interference of tree and two-loop graphs can 
be reconstructed. In previous work [1, 2], the finite part has been separated from infrared singular 
parts using the Catani formalism [49]. In this procedure, the singular operators multiply tree and 
one-loop amplitudes. Changing the scheme for the external particles therefore changes the overall 
pole contribution. However, the finite remainder left after subtracting the poles is invariant under 
changing from the conventional dimensional regularisation scheme to the 't Hooft-Veltman scheme. 
Therefore, the finite parts of the helicity amplitudes presented here precisely reproduce the finite parts 
of the interference of tree and two-loop graphs given in [1, 2]. Similarly, the finite parts of the self 
interference of one-loop amplitudes can be reconstructed purely from the finite parts of the one-loop 
helicity amplitudes. In addition, the finite remainders of the two-loop amplitudes for quark-gluon 
scattering in the HV scheme agree with the independent calculation of Ref. [3] . 

For the case of quark-gluon scattering, the presence of four coloured external particles gives rise 
to additional colour correlations proportional to 1/e. While the existence of such additional non- 
trivial colour correlations was pointed out in Ref. [49] and are expected on general grounds [57], the 
precise form of the colour structure of H^^ 2 ' was not predicted. So far, the results for gluon-gluon 
scattering [30] as well as those for quark-gluon scattering presented in this paper and in Ref. [3] , support 
the form given in Eqs. (3.38) and (3.39). These colour correlations vanish when the interference of 
tree and two-loop amplitudes is constructed. 

The more general goal is to evaluate the NNLO QCD corrections to the associated scattering 
processes which should yield a significant improvement in theoretical stability, see for example [4]. 
However, the reality is that the phenomenology is still some way off and there is still a significant 
amount of work remaining to be done. The two-loop corrections to the 2^2 partonic process 
presented here, need to be combined with infrared singular contributions from both the one-loop 2^3 
partonic process, when one of the partons becomes unresolved [65, 66, 63] and the tree-level 2^4 
partonic process when two partons become unresolved [67, 68, 69, 70, 71] . These different contributions 
need to be put together in such a way to make the explicit singularities analytically cancel, so that 
the four-dimensional limit can be safely taken, and then by encoding the finite remainders in a stable 
numerical implementation. There has been recent progress in identifying appropriate counterterms [72, 
73] and we are optimistic that the technical problems will soon be overcome enabling predictions of 
both dijet, jet and photon and diphoton production at the Tevatron and the LHC at an accuracy of 
0(5%) which is expected to be competitive with the experimental measurements. 
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A. Finite two-loop contributions 

In this appendix we give explicit formulae for the coefficients A,...,L for the qq — > gg process. 
Expressions for quark-gluon scattering as well as the processes involving photons are available from 
the authors as a FORM output file. 

As usual, the polylogarithms Li n (w) are defined by 



f w dt 

Li„(u>) = / — Li n _i(i) for n = 2,3,4 
Jo t 
t' w dt 

U 2 (w) = - -log(l-t). (A.l) 
Jo 1 



Using the standard polylogarithm identities [74], we retain the polylogarithms with arguments x, 1 — x 
and (x — l)/x, where 

t u u x — 1 

x = — , y= = l-x, z = -- = . (A.2) 

s s t x 

For convenience, we also introduce the following logarithms 

* = 10g(^), >' = lo«(v). * = 

where \i is the renormalisation scale. 

A.l q(p 2 , +) + q(pi, -) -> g(p 3 , +) + g(p4, +) 
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+ 1 X 4 - 2 X 5 Y + - X 3 + X 2 Y 2 - - X 2 Y + 1 X 2 ir 2 + - X 2 - -XY 2 

6 3 12 12 3 9 6 

2 9 7 1 , 1 9 7 9 313 

7 4 85 2 4085\ 
+ 45" "72" + 1296 J 

+ f - 12Li 4 (x) + 12Li 4 (y) + 12Li 4 (z) + 12Li 3 (x) F + 8Li 3 (y) - 2Li 2 (x) vr 2 

-8 Li 2 (y) y + g ^ 4 - 2 A 3 K + 1 ^ 3 + 3 X 2 y 2 - - X 2 Y + X 2 vr 2 + 2 X 2 - 4 X Y 2 
9 2 5 4 9 7 4 7 9 \ i 

_2xy7r 2 --xy + -x-i2yc3 + ^y^ 2 -8C 3 + 77^-0^ - 

3 3 3 15 9 / u 

+ ^ - 8 Li 4 (x) + 8 Li 4 (y) + 8 Li 4 (z) + 8 Li 3 (x) V - Li 3 (x) + ? Li 2 (x) X 

Au 2 (x)k 2 + - X* -- X 3 Y + 2X 2 Y 2 + - X 2 k 2 + - X 2 -- XY v 2 + - X<k 2 

3 w 3 3 3 23 9 

49 9 9 2 \ s 

+ 12 3^ 3 )^ 



3 ,, 2 s 



2 



+ 2 y * 

+ ^4Li 3 (x) + |Li 2 (x) + ^ 2 + ^y + ^-2y-4C3 + ^ 

+ ^12Li 3 (x) + 8Li 2 (x) + ^ 2 -^y+^X-^y-12C3 + 3^ * 
+ ( 8Li 3 (x) + 2 - Li 2 (x) - ^ 2 + X - 8C 3 - I vr 2 j ^ 
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(A.34) 

F^M = ( —Y + —} iir + ( - - Y 2 + — Y + — 7T 2 } 
b+ - + - yi8 27 I \ 12 27 27 J 

(A.35) 

« - (-l*> + l*r-ly») + (Ix-l^.f^lxr-l^) 1 

(A.36) 



- 4 Li 4 (x) - - Li 4 (y) + 4 Li 4 (z) + Li 3 (x) X + - Li 3 (x) Y - — Li 3 (x) - - Li 3 (y) X 

+ y Li 3 (y) y + % - Lia(y) - Li 2 (x) X 2 + y Li 2 (x) X + y Li 2 (x) vr 2 + 2 Li 2 (y) X Y 

- ? Li 2 (y) y 2 - ? Li 2 (y) y + ^X 4 - X 3 Y + 1 X 2 y 2 + \x 2 Y - -X 2 vr 2 - \ X Y 3 
4 4 o 4 3 6 4 

23 9 11 5 91 , 1 4 11 o 1 o 9 161 9 

-^Ay 2 + -xy + -xc3-^x. 2 --y^ + -y3- s y 2 7r 2 + -y 2 -3y C 3 

43 2 211 v 9 A 43 4 391 2 \ 
+ 48^ -108 y -4 C3 + 120 7r " re-^J 

— Li 4 (x) - - Lia(x) X - Lia(x) K - y Li 3 (x) - 2 Li 3 (y) X - - Li 3 (y) 
+ ^ Li 2 (x) X 2 + y Li 2 (x) X + | Li 2 (x) vr 2 + Li 2 (y) IF + y Li 2 (y) y + ^ X 4 

+^x 3 y-^x 3 + |x 2 y-^x 2 vr 2 + ^x 2 + ^xy 2 + ^xy7r 2 + ^xy 

6 12 6 8 24 4 3 4 

+ 2^C3- ¥ Xvr -_X + y < 3__y 7r + y C 3 - ^ vr - - tt j - 

( 33 33 7 31 31 

11 Li 4 (x) - — Li 4 (y) - y Li 4 (z) + - Li 3 (x) X-— Li 3 (x) Y + — Li 3 (x) 

o 91 in 11 

+2 Li 3 (y) X - - Li 2 (x) X 2 - - Li 2 (x) X + - Li 2 (x) vr 2 - Li 2 (y) X Y - - X 4 
4 12 11 lb 

+§x 3 y-lx 3 -f x 2 y 2 + ix 2 y-Ix 2 vr 2 + ^x 2 + |xy 7 r 2 

7 ^ 191 r^, 31 31 139 A i 2 

"S^ 6 " 72** + Y yC3 -12 C3 - 360 7r 

(55 13 31 

- Li 4 (x) - - Li 4 (y) + 8 LU(z) - y Li 3 (x) X + 6Li 3 (x) y - — Li 3 (x) 

13 31 5 31 11 

-6Li 3 (y) X + — Li 3 (y) y + — Li 3 (y) + - Li 2 (x) X 2 + — Li 2 (x) X + — Li 2 (x) tt 2 
l b 4 b 3 
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+4 U 2 (y) XY-\ U 2 (y) Y 2 - f Li 2 (y) Y + f g X 4 - ^ X 3 Y - H X 3 + 3 X 2 y 2 

7 o 5 9 n 43 o 1 , 7 n 5 „ 1 11 

4 12 18 12 4 6 2 2 S3 

23 n 19 1 A 11 o 1 9 9 26 9 11 23 o 

"16 X * + 27 * " 16 5 + 18 F + 4 Y * ~ IT 1 " T F 6 + IB y ' 
_19 «_ 1 \t 3 s 
27 180 A j s 2 t 

+ Li 3 (z) + 5 Li 3 (y) - 4 Li 2 (z) X + g Li 2 (x) - | Li 2 (y) F + * X 3 - * X 2 Y 

1 o 19 1 oH 1 o 19 n 3 n 149 1 

+ A X + ^ XY ~2 X * + Y X ~4 Y ~2A Y ~2 Y * + 12 Y ~2^ 
_187 2 _21l\ 
144 108 j 

(g -j g-i w frr 

- - Lia(x) - 2 Lig(y) - Li 2 (x) X - — Li 2 (x) + Li 2 (y) F + - X 3 - X 2 + - X Y 

1.2 19, 17 y , 9 , 5 2 5\ t 
"3 X7T + Y X + T y+ 2 C3 -9 7r "ej n 

+ ^ - 12 Li 3 (x) + 2 Li 3 (y) - \ Li 2 (x) X - ^ Li 2 (x) - Li 2 (y) y - I X 3 + ^ X 2 
+ ( - ^ Lig(x) + \ Li 3 (y) - | Li 2 (x) X + | Li 2 (x) + ? Li 2 (y) Y + ^ X 3 

- 5 x 2 y- 5 x 2 + 5 xy 2 + 31 xy + xvr 2 - 5 x--y 3 + 5 y 2 -y7r 2 

4 6 4 6 6 12 6 

5., 31 ,\ i , A 

+ l Y -m")- s +3Y t)"' 

(A.37) 

= ( - 2 Li 4(*) " f LUG/) + Lia(x) X + 1 Li 3 (z) y - | Li 3 (x) + 1 Li 3 (y) X + | Li 3 (y) Y 
+ 9 - Li 3 (y) + | Li 2 (x) X - I Li 2 (s) vr 2 - 1 Li 2 (y) Y 2 - \ Li 2 (y) y - ^ X 3 Y 

+^x 2 y 2 -^x 2 y + ^x 2 7r 2 -^xy 3 -^xy 2 -^xyvr 2 + ^xy 

1 4 3 2 \ 
36 8 i 
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(23 11 
- — Li 4 (x) - 3 LU(y) - 3 U 4 (z) + — Li 3 (x) X - Li 3 (x) F - 4 Li 3 (x) + 4 Li 3 (y) X 

+ * Li 3 (y) - | Li 2 (x) X 2 + 4Li 2 (x) X - ^ Li 2 (x) vr 2 - 2 Li 2 (y) 17-^ Li 2 (y) F 

_A x 4 + -x 3 F--x 2 F 2 --x 2 F + x 2 7r 2 + -x 2 --XF 2 --XF7r 2 

32 6 4 12 8 4 6 

12 2 Sd 3 2 S 12 360 24 J u 

+ ^ - 4Li 4 (x) - | Li 4 (y) - ^ LU(z) + ^ Li 3 (x) X- 1 - Li 3 (x) F - ? Li 3 (x) + 2 Li 3 (y) X 

Li 2 (x) X 2 + ? Li 2 (x) X - ^ Li 2 (x) vr 2 - Li 2 (y) XF-Ax 4 + lx 3 F-^X 3 

_|x 2 F 2 + -X 2 F + -X 2 vr 2 + Ix 2 --XFvr 2 - -X Cs + - X vr 2 
8 3 8 8 12 2 S3 8 

+ 2 yC3 + 4 C3 + 360 7r J^ 

+ ^ Li 4 (x) - l - Li 4 (y) - 1 Li 3 (x) X+ l - Li 3 (y) F + 1 Li 2 (x) X 2 - 1 Li 2 (x) vr 2 
-^ i2(y )F 2 -l^ + lx3F + lx 2 . 2 + ix 2 -lxF3-IxF. 2 

+ ^-f^ + ^ 4 -l^^ 2 ^-^ 2 + f^^ 4 -^ 2 )^ 

l 9 s/l 9 l lolo\t 2 

4 i I 4 2 4 4/s 2 

+ ^ (j Li 3 (x) + 2 Li 3 (y) + ^ Li 2 (x) - 1 Li 2 (y) Y - \ X 3 + \ X 2 Y - \ X 2 

+ (j Li 3 (x) + 4 Li 3 (y) - I Li 2 (x) X + ? Li 2 (x) - 2 Li 2 (y) F - | X 3 - ^ X 2 - 1 X F 

1 2 8 7 9 1 2 l\i 
-12 X7r + 3 X -12 y -2 C3 -2^ + 2"k 



+ ( :; Li 3 (x) + 2 Li 3 (y) - l - Li 2 (x) X + *j Li 2 (x) - U 2 (y) Y - 1 X 3 - || X 2 + | X F 



2 



+ ( - \ Li 3 (x) + 1 Li 3 (y) + 1 Li 2 (x) X - 1 Li 2 (y) F - 1 X 3 + 1 X 2 F 
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XY 2 + 2X + ^-Y 3 - 2 y J - + l ) Y-)iTT 



(A.38) 



J!f"l' m = Lia(x) + ^ Li3(y) - \ Li 2 (x) X - ± Li 2 (y) Y - A X 2 Y + 1 X y 2 

29 5 9 1 , 13 n 17 n 19 1 17 o\ 

— xy + — Xvr 2 --y 3 — y 2 y^ 2 + — y- -Cs + — vr 2 

72 12 9 24 72 72 3 s 12 J 

+ (- LU(x) + Li 4 (y) + Li 4 (z) + Li 3 (x) y + 2 Lig(y) - ^ Li 2 (x) vr 2 - 2 Li 2 (y) Y 

+ Y4 XA -1 X3Y+ 1 X3 + 4 X2Y2 -T4 X2Y+ ^ X27t2 -I2 X2 - XY2 
1 o 1 3 9 1 1 9 7 4 5 9 \ i 

--xy 7r 2 + -xy + -x 7 r 2 + -x-y C 3 + 3y. 2 -2 C 3 + ^. 4 -^ 7 r 2 j- 

+ ^ - 2 Li 4 (x) + 2 Li 4 (y) + 2 Li 4 (z) + 2 Li 3 (x) Y - ^ Li 2 (x) vr 2 + ^ X 4 - ^ X 3 y 
+ 1 x 2 Y 2 + i- X 2 y + - X 2 vr 2 + \ X 2 - - X Y vr 2 - — X vr 2 - 2 y Cs + - vr 4 > 



2 24 6 4 3 12 ' 90 I u 

+ ^ Li 3 (x) - I Li 3 (y) - I Li 2 (x) X + \ Li 2 (y) y + 1 X 3 - 1 X 2 Y - | X 2 

+-xy 2 + xy + -X7r 2 + -x--y 3 + -y 2 --y7r 2 --y--7r 2 ^- 

4 72 18 9 24 72 18 2 js 

2 t 



2 



1 „o 1 1 ,,o 1 2 \ i 



+ --x z + -xy--y z --7r 



2 



S 2 



+ -iLi 2 (x ) + ix 2 -^y-^x + A y2 _I! y + iZ 7r2 + ^ 

11 3 w 8 3 72 24 24 72 72 

+ + 2 Li 2 (x) - 1 a" - 1 x r - I A- + 1 y - c, + i J 1 

I 3 v 7 12 3 36 12 36 18 / s 

-Y^Ur 



(A.39) 
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K g 9 (2)[3] = I _ 1 ( ) + I Li / )Y+ ±X 2 Y + ^XY 2 +-XY--Xk 2 -- Y 2 +- Y tt 2 
s+-+- I 3 -i\yj t 3 i\V) -r ^ -r 24 -r g 12 2 36 

+ 24 y+ 3 C3 "8 7r J 

+ ^ - 2 Li 4 (x) + 2 Li 4 (y) + 2 Li 4 (z) + 2 Li 3 (x) F + X - Li 3 (x) + 4 Li 3 (y) - X - Li 2 (x) X 

- ^Li 2 (x) vr 2 - 4 Ii 2 (y) y + 1 X 4 - 1 X 3 y + ^x 2 y 2 + ^ * 2 Y + 1 X 2 vr 2 

+ -X 2 -2Xy 2 --Xy7r 2 --Xy--Xvr 2 + -X-2yC3 + -Fvr 2 
4 3 4 9 6 s 3 

13 ^ + 7 ^ 17 ,\ i 

3 90 36 I u 

(11 2 
- 4Li 4 (x) + 4Li 4 (y) + 4Li 4 (z) + 4Li 3 (x) F - - Li 3 (x) + - Li 2 (x) X - - Li 2 (x) vr 2 

+ -X 4 - -X 3 Y + X 2 Y 2 + ^rX 2 y + -X 2 vr 2 + \x 2 - -XYir 2 + — X tt 2 
6 3 24 3 4 3 36 

+ ^ Li 3 (x) - 1 Li 3 (y) - 1 Li 2 (x) X + l -U 2 {y)Y +\x 2 - l -XY - ^Xir 2 - l -X 

1 n 5 o 1 1 ,\t 

+ i y + i8 y7r + 3 y+ r J« 

4 t 

+ ff_iu j(l) + ijr»_i jr y + ijr-iy»-2y + » , \ 

1 I 3 w 24 6 8 8 8 24 J 

+ (^i3(x) + iLi 2 (x)-I^ + ixy + ix-4 C ,-i^^ 

3 V "V 
Y - \ ITT 

2 t I 



(A.40) 



L^i 2 i'L 3] = I -IxY-Iy 2 + ^Y-Ik 2 ) +(Ix 2 -^X-Iy 2 + Iy X 1 



72 36 54 72 / V 36 54 36 54 / s 
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18 18 J s j 

(A.41) 

A.3 q(p 2 , +) + q(pi, -) -► $(p 3 , -) + ff(P4, +) 

The finite remainders are obtained by Bose symmetry from the g(j>2, +)+<j(pi, —) —> g(P3, +) + <?(P4 ) — ) 
process. 

Af/™ =A 9 s ^3(t~u), (A.42) 

etc. 

A.4 q(p 2 , +) + q(pi, -) -> #(p 3 , -) + ff(p4, -) 

The finite remainders are obtained by Bose symmetry from the q(p2, +)+<j(pi, —) —> g(P3, +) + <?(P4 ) — ) 
process. 

= A 9jmm {t ^ u) ^ (A .43) 

etc. 



- 41 - 



B. Finite one-loop contributions 



In this appendix we list the coefficients A, . . . , E for the finite one-loop amplitudes defined in Section 4.3 
for qq — > gg. Expressions for quark-gluon scattering as well as the processes involving photons are 
available from the authors as a FORM output file. 

B.l q(p 2 , +) + q(pi, -) -> g(p 3 , +) + g(p4, +) 



4 <?<?,(l),[l] 




1 s 

2 t 


1 

3 


Dflfl.(l).[l] 

-°s+-++ 




1 s 

2 t 






1 

~ 3 






4flfl,(l),[2] 

^3 + - + + 




1 t 

2 u 


1 

6 


Dflfl,(l),[2] 




1 t 

2 u 


1 

" 2 






1 

3 






= 







B.2 g(p 2 , +) + ?(pi, -) -> 5(P3, +) + S(P4, -) 



3 n 


5 






2* + 




x + 


- - 




2 


2 y « 


5\ t 








2"h + 


( 


3X 













• :! ) +f3^ + |)-+f3X + l)^ + X^)^ (B.3) 



C^ ] = -^X+iiTT (B.l) 



5' a i'I 2l = [ -^ + ^^ + 3) -\y^-+{ { -2Y-1) -y S )zvr (B.5) 



3 / 2 t 1 1 3 t 



B 99,(1)P] = _ A __ Y ^__ Ym ~_ (B6) 



1 



cfi ,(l),[2] = _ - y - 



1 „ 1 



3 3 



(B.7) 



D 99,(i),m = f xy _i y2 + o y _i 7r2 + 3 X 



2 



's+-+- 



2/4 « 
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2 u 



5 1 „ 5 
X - - Y 2 + - Y 



+\x-+(x-Y\-\m (B.8) 





B.3 q(j)2, +) + q(pi, -) -> #(p 3 , -) + g(p4, +) 

The finite remainders are obtained by Bose symmetry from the q(p2, +)+q(pi, —) —> g(P3, +) + <?(P4 ) — ) 
process. 



^sH h ~~ A s-\ 1 — K L 



AfMM = A f^l [ 3 (<<-«), (B.10) 



etc. 



B.4 q(p 2 , +) + q(pi, -) -> 5-(P3, -) + s(p4, -) 

The finite remainders are obtained by Bose symmetry from the q(p2, +)+q(pi, —) —* g(P3, +) + <?(P4 ) — ) 
process. 

A 9£WM = Af^Wf (/ <- u), (B.ll) 

etc. 
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